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Abstract 

We consider the dynamic escape of a small packet of ultracold atoms launched from within an 
optical dipole trap. Based on a theoretical analysis of the underlying nonlinear dynamics, we predict 
that fractal behavior can be seen in the escape data. This data would be collected by measuring 
the time-dependent escape rate for packets launched over a range of angles. This fractal pattern 
is particularly well resolved below the Bose-Einstein transition temperature — a direct result of the 
extreme phase space localization of the condensate. We predict that several self-similar layers of 
this novel fractal should be measurable and we explain how this fractal pattern can be predicted 
and analyzed with recently developed techniques in symbolic dynamics. 

PACS numbers: 32.80.Pj, 05.45. Ac, 05.45.Df 
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Chaotic escape is a widespread transport process underlying such diverse phenomena as 
conductance through ballistic microstructures 1|, emission from deformed micro-disk semi- 
conductor lasers 2L molecular scattering and dissociation {jj, celestial transport \&\, and 
atomic ionization 3, (|. We have been particularly motivated by the chaotic ionization 
of hydrogen in applied parallel electric and magnetic fields, for which a recent theoretical 
analysis [6j] predicts that the time-spectrum for ionization will display a chaos-induced train 
of electron pulses. This prediction is based on classical ionizing trajectories (Fig. [Ik^, which 
propagate from the nucleus into the ionization channel, via the Stark saddle. These tra- 
jectories exhibit fractal self-similarity, which is reflected in the pulse train. This prediction 
has been recently confirmed by full quantum computations 7]. However, the experimental 
observation of these chaos-induced pulse trains remains unrealized. 

In this Letter, we propose an alternate physical system — the escape of ultracold atoms 
in specially tailored optical dipole traps — that exhibits similar escape dynamics. However, 
we show that the flexibility and control afforded by cold atoms, especially in engineering 
the initial state, should permit the direct imaging of fractals in the escape dynamics, in- 
cluding novel self-similar features. Furthermore, we show how this self-similarity can be 
analyzed using a recently developed symbolic formalism. Fundamentally, these fractals re- 
sult from homoclinic tangles [8y — a general mechanism for phase space transport. Hence this 
Letter suggests that cold atoms could serve as a unique high-precision experimental probe 
of this mechanism. Finally, the cold atom experiments discussed here are readily feasible 
with present-day experimental configurations and should prove easier to realize than the 
previously mentioned ionization experiments. 

Recent experiments from the Raizen 9] and Davidson lOj groups have made first steps 
along these lines. They independently measured the long-time survival probability for ul- 
tracold atoms escaping through a hole in an optical billiard, demonstrating the distinction 
between regular and chaotic escape dynamics. In contrast, our Letter focuses on the short- 
to intermediate-time dynamics, where fundamentally distinct phenomena, such as fractal 
self-similarity, are predicted to appear. 

The double- Gaussian trap: We consider a dipole potential consisting of two overlapping 
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Gaussian wells 



V(x,y) = -V 1 exp(-[(x/ ( T lx .) 2 + (i/M,) 2 ]/2) 

- V 2 exp{-[({x - x 2 )/a 2x f + (y/a 2y ) 2 }/2), (1) 

as shown in Fig. [lb. This potential can be created by two red-detuned, far-off-resonant 
Gaussian beams; atomic motion can be further restricted transverse to the xy-pl&ne by 
a uniform laser sheet. Here, we take a± x = 0.18, a± y = 0.24, a 2x = 1-08, o 2y = 0.24, 
x 2 = 0.72 (measured in millimeters), and V\ — V 2 — 35.5 (measured in recoil energies 
E T = h 2 kl/2m Kh = h ■ 3.77 kHz, where A = 2Tr/k L = 780.2 nm for the D 2 transition of 
87 Rb.) The double-Gaussian potential shares several features in common with the hydrogen 
potential in Fig. [Th- The "primary" Gaussian centered at the origin is analogous to the 
Coulomb well; the elongated "secondary" Gaussian on the right is analogous to the ionization 
channel; and the saddle connecting the two Gaussian wells is analogous to the Stark saddle. 
We are interested in the transport of atoms from the primary well into the secondary well. 
Fig. [lb shows two representative trajectories that move away from the origin with initial 
speed 4.12 cm/s, pass over an unstable periodic orbit (PO) near the saddle, and then strike a 
resonant laser sheet (the vertical dashed line.) This sheet forms a detection line that serves 
both to image the escaping atoms and to scatter them out of the trap, preventing their 
return into the primary well. 

Figured plots the time for a trajectory with initial speed 4.12 cm/s (energy -14.9 D 2 
recoils) to move from the origin to the detection line, as a function of launch angle 9, 
measured relative to the positive x axis. The resulting escape-time plot is highly singular, 
with numerous icicle-shaped regions whose edges tend toward infinity. These "icicles" exhibit 
a self-similar fractal pattern. Such patterns also occur in the chaotic ionization of hydrogen 
and are characteristic of chaotic escape and scattering. 

A proposed experiment to measure self-similar patterns in the escape-time plot: We con- 
sider a small Gaussian packet of ultracold atoms launched from the origin with speed 4.12 
cm/s and 9 = 2.04 (the right line in Fig. [2h). The subsequent flux of atomic trajectories at 
the detection line is then computed as a function of time. Fig. [3h shows this flux for an initial 
thermal packet that occupies a phase space area 500 times Planck's constant in both the x 
and y degrees of freedom. Fig. [3b, on the other hand, uses a packet that occupies a single 
Planck cell, appropriate for a pure dilute Bose-Einstein condensate (BEC) in the regime 
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of negligible interactions. The condensate packet closely follows the trajectory in Fig. [Tb, 
exiting as a sharp pulse at 160 ms, near the bottom of the 9 = 2.04 icicle in Fig. [2^. The 
thermal packet also produces a pulse at 160 ms, but its larger phase space extent populates 
neighboring icicles, thereby producing the additional pulses in Fig. [3^. 

By repeating the preceding computation for different launch angles of the thermal packet, 
we obtain the aggregate data in Fig. [2b, where the shading records the atomic flux as a 
function of arrival time and the packet's launch angle. Fig. [3ji then corresponds to the 
vertical slice through Fig. [2b at 9 = 2.04. The thermal data appear as a blurred version 
of the sharp escape-time plot. For example, the left and right icicles are associated with 
prominent dark patches, and in between a few wispy patches can be associated with the 
bottoms of other icicles. Overall, however, the intricate icicle structure is poorly resolved 
by the thermal data. 

A remarkable increase in resolution occurs below the BEC transition, shown in Fig. H}. 
Many icicles are now clearly resolved, a direct consequence of the high phase-space localiza- 
tion afforded by the condensate. This increase in resolution prompts us to look deeper into 
the fractal, with the expectation that we can directly measure its fine-scale structure. 

Self- similarity of the escape-time data: The first column of Fig. H] shows data plotted for 
three distinct intervals of the launch angle. (We concentrate at present on the upper plot in 
each pair.) The three plots look remarkably similar, and icicles in one plot can be identified 
with icicles in the other two. In fact, this pattern of icicles occurs throughout the escape- 
time plot and on all scales. (See below.) One of the principal observations of this Letter is 
that these structures are resolved by the overlaid BEC data. The icicles look progressively 
more blurred as we move down the column because the interval width is decreasing and the 
escape time is increasing. 

The pattern in the first column of Fig. H] is not the only repeated pattern. The second 
column of Fig. [4] shows another pattern that also occurs on all scales throughout the same 
escape-time plot. As we will see, many such repeated patterns, or templates, exist in the 
escape-time plot. Within a given template, all other templates can be found on smaller 
scales. That is, each template occurs as a subtemplate of every other template in an infinitely 
recursive nesting. Our computations predict that several nested layers will be experimentally 
visible. 

Experimentally, the observation of these phenomena will not be easy but certainly fea- 
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sible. With a 1.06 jim. laser, the above trap geometry can be realized with about 80 W 
of power. With this detuning, a 87 Rb atom has at most a 4% probability of spontaneous 
scattering over a half second. Acceleration of the atoms to the initial velocity of 4.12 cm/s is 
easily accomplished; for example, a chirped, one-dimensional optical lattice of 785 nm light 
with 50 mW of single-beam power and beam radius Wq = 100 /im can accelerate the atoms 
in about 1.6 ms with negligible (< 1%) probability of spontaneous scattering and heating 
due to energy-band transitions. The primary difficulty is the subrecoil initial conditions 
required, even in the thermal case. However, a standard expanded BEC should suffice. 

Theoretical foundations of the fractal structure: We next describe how the self-similar 
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fractal data can be described with recently developed symbolic tools [11]. We first specify a 
two-dimensional surface of section in the xypxPy-pha.se space by fixing the energy E at -14.9 
D 2 recoils and setting y — 0. Thus, every time a trajectory passes through the x-axis, we 
can record (x,p x ), defining a Poincare map that maps a given intersection (x,p x ) forward 
to the next intersection (x',p' x ). Fig. [5] shows the corresponding surface-of-section plot. The 
vertical line Co at x = consists of all points at the origin moving outward with arbitrary 
launch angle 9. The point zx is the unstable fixed point equal to the intersection of the 
unstable periodic orbit (PO) in Fig. [Tb with the surface of section. Attached to zx are its 
stable S (thick) and unstable U (thin) manifolds, consisting of all points that asymptote to 
zx in the forward and backward directions, respectively. These manifolds intersect an infinite 
number of times, forming an intricate pattern called a homoclinic tangle flfl Q,Q- The 
segments of the stable and unstable manifolds connecting zx to the point Po define the 
shaded region called the "complex". 

Escape from the complex occurs via escape lobes E n , defined in Fig. [5] as the regions 
bounded by the stable and unstable segments connecting P n to Q n . The lobe E_i, inside 
the complex, maps to E , outside the complex. Once in E a point then maps to E%, E 2 , etc., 
eventually passing into the secondary well on the right. The lobes contain all points 
that escape in k iterates. These lobes become progressively more stretched and folded as 
k increases. (An analogous sequence of lobes C n controls capture into the complex.) Note 
that we are able to chose physical parameters that make the lobes quite large compared to 
Planck's constant (Fig. EJ). 

We plot as a function of 9 the number of iterates for a point to escape the complex, shown 
as the lower plot of each pair in Fig. HI These plots straighten each icicle into a constant 
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escape segment. A segment that escapes on iterate k is an intersection between and Cq. 



For examp 
Ref. 



e, the segment at iterate two in Fig. H} is the intersection with E_ 2 in Fig. [5j 



11] introduces a symbolic technique, called homotopic lobe dynamics, to compute 
the structure of escape segments based on the tangle topology. (See also Refs. [jjj].) We 
summarize the results obtained from applying this technique to the tangle in Fig. 

The structure of the discrete-escape-time plot up to a given iterate n is specified by a 
string £ n of symbols in the set {ci, c 2 , a, b, f, w 0) v>i, u 2 , • ••} as well as their inverses, e.g. of. 
The first string is l\ = Mi6 _1 m . All subsequent strings can be obtained from the first by 
mapping each symbol forward according to the substitution rules: 

ci c 2 , c 2 ^ f' l u Q auQ V> (2a) 

a 1- ft" V&> / ^ cr^o 7, (2b) 

Mn^«n+l, b^b^UQ 1 /, (2c) 

using the standard convention for iterating inverses, e.g. 1— > f~ l Uob. For example, the 
first four strings are: 

i\ = wi& _1 Mo, (3a) 

h = u 2 f" l Uobu 1 , (3b) 

h = u 3 f' 1 uoc 1 u 1 b' 1 UQ 1 fu 2 , (3c) 

£ 4 = U4f~ 1 ^c 1 u 1 c 2 u 2 f' 1 uobui 1 Ci 1 u^_fu 3 . (3d) 

An £ n string encodes the discrete-escape-time plot as follows. Each appearance of u^ 1 in £ n 
(underlined for emphasis) represents a segment that escapes at iterate n. For example, the 
u Q factor in Eq. (I3"a|) corresponds to the escape segment at iterate one in Fig. Hi. This u 
factor then maps forward to u% in Eq. (I3bj) . In general, we see that each u^ 1 in £ n represents 
a segment that escapes at iterate n — k. In Eq. (pbj) . another m factor has appeared, 
corresponding to the escape segment at iterate two. This segment is to the left of the first 
segment, just as the uq factor in Eq. fl3b|) is to the left of the U\ factor. In general, the 



left-right ordering of 1 symbols in £ n represents the left-right ordering of segments in the 
discrete-escape-time plot. On the next two iterates, two new u factors appear in Eq. ([5c]) . 
and three more appear in Eq. (13dj) . in agreement with Fig. Hi. 

All other symbols besides u^ 1 represent gaps between adjacent escape segments. For 
example, b~ l in Eq. fl3al) represents the gap b~ l {i) in Fig. H}, and b~ l in Eq. fl3c|) represents 



the gap b~ l (ii) in Figs. H} and Hfr. The string £5 will also contain a 6 _1 factor, representing 
the gap h~ l (iii) in Figs. H^i andH}ii. Since each 6 _1 factor generates exactly the same string 
of symbols under Eqs. each gap labeled b~ l in Fig. @] contains the same pattern of escape 
segments. This means that b^ 1 corresponds to a particular template, i.e. to find occurrences 
of this template in the escape-time plot, we need only look for occurrences of b~ l in the 
expression for £ n . 

It follows from the above logic that each symbol {ci, C2, a, b, /} generates its own template, 
with inverse symbols generating reflected templates. For example, mapping / _1 forward 
three times yields 

f^Uo^U^Uif^UoCLUQ 1 /. (4) 

The reader may verify that Eq. (jlj) describes the segments up to iterates five, seven, and 
six in Figs. H^i, Hb, and Hfc. Note that different experimental parameters will yield different 
algebraic rules and different templates. 

This algebraic formalism computes a minimal set of escape segments, but generally not 
all segments. That is, at later times, we typically find additional segments in the numerics. 



This illustrates what has previously been called an epistrophic fractal 13j. Nevertheless, 
unpredicted segments can be accommodated within an updated algebraic formalism, as 



explained in Ref. 



Ill- 



Conclusions: We predict that experiments on the intermediate-time escape dynamics of 
ultracold atoms from an optical trap can directly image fractals. The resolution is particu- 
larly good when using a BEC. The fractal structure depends on the topology of homoclinic 
tangles, which are common to numerous chaotic systems. Such experiments would thus pro- 
vide a new laboratory tool for the study of an important chaotic mechanism. Similarly, an 
improved understanding of the chaotic escape pathways of atoms from optical traps could 
be relevant for the understanding of mixing and thermalization in traps and for the control 
and coherent emission of atomic wavepackets. Finally, the dependence of these fractals on 
atom density could serve as an interesting probe of atom-atom interactions, a subject to be 
explored in future work. 
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FIG. 1: a) The potential energy for a hydrogenic electron in applied parallel electric and magnetic 
fields. Two ionizing trajectories are shown, b) The double-Gaussian potential. Gaussian widths 
are indicated by the dashed ellipses. 
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FIG. 2: a) The escape-time plot. Vertical lines denote 9 for the trajectories in Fig. b) The 
shading shows the escaping flux from an initial packet of size 500 times Planck's constant. 
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FIG. 3: The flux of atoms reaching the detection line for two different Gaussian ensembles launched 



at 9 = 2.04. The position and momentum widths of the initial ensembles are given in units of mm 
and recoil momenta (p T = hk L = m Rb • 5.88 mm/s for the Rb D2 transition.) 
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FIG. 4: Escape time over six different angular intervals. Each escape-time plot is matched below 
by its corresponding discrete-escape-time plot. 
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FIG. 5: The surface of section, showing the unstable fixed point zx and its associated homoclinic 
tangle. The left rectangle has area ten times Planck's constant. 
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